Correlations and thermalization in driven cavity arrays 
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Abstract. We show that long-distance steady-state quantum correlations (entanglement) between 
pairs of cavity-atom systems in an array of lossy and driven coupled resonators can be established 
and controlled. The maximal of entanglement for any pair is achieved when their corresponding 
direct coupling is much smaller than their individual couplings to the third party. This effect is 
reminiscent of the coherent trapping of the A— type three-level atoms using two classical coherent 
fields. Different geometries for coherent control are considered. For finite temperature, the steady 
state of the coupled lossy atom-cavity arrays with driving fields is in general not a thermal state. 
Using an appropriate distance measure for quantum states, we find that the change rate of the 
degree of thermalization with respect to the driving strength is consistent with the entanglement of 
the system. 

PACS numbers: 03.67.Bg, 03.67.Hk, 03.65.Yz, 42.50.Pq 



I. INTRODUCTION 

Coupled cavity arrays have recently been proposed as a novel system for realizing quantum computation [l| and 
for simulations of quantum many-body systems 0. More recently, the steady-state polaritonic Q and membrane 
entanglement [2] of driven cavity arrays were studied under realistic dissipation environment. Also, there has been an 
attempt to relate coupled cavity arrays with Josephson oscillations ■ 

At finite temperature, it is expected that the steady state of the coupled-cavity system is a thermal state, since 
standard statistical mechanics tells us that if a system interacts with a large reservoir at a fixed temperature, it will 
relax eventually to an equilibrium state characterized by the Boltzmann distribution with a well-defined temperature, 
i.e. that of the reservoir. However jSuch thermal relaxation is true only for some simple systems such as a single empty 
cavity coupled to a thermal bath [6|. For many other systems e.g. coupled cavities with external pumping lasers, the 
steady state does not need to be a thermal state, and its deviation from a thermal state depends on various factors: 
inter-cavity couplings, presence of the pump, detuning and so forth. 

The purpose of this article is twofold: firstly we wish to demonstrate the possibility of achieving coherent control of 
the steady-state entanglement between mixed light-matter excitations generated in macroscopically separated atom- 
cavity systems, and secondly, we hope to elucidate the conditions under which the steady state differs from a thermal 
state, especially the relation between the thermalization of the system and the correlations of the subsystems, using 
the coupled atom-cavity system as an example. This paper is organized as follows. In Sec. II, we introduce the setup 
and the Hamiltonian for coherent control of the steady-state entanglement. In Sec. Ill, we derive an effective equation 
for the dynamics of the system. In Sec. IV, we discuss the coherent control of the steady-state entanglement. In 
Sec. V, we discuss an alternative setup: two coupled cavities with three driving fields. In Sec. VI, we discuss the 
thermalization of two defect cavities coupled to one driven wave guide in between. In Sec. VII, we summarize our 
result. 



II. THE SETUP AND THE HAMILTONIAN 



The setup we study is shown in Fig. [T] It contains three interacting atom-cavity systems {Si, S2, S^) connected 
by three waveguides/fibres. Each waveguide/fibre is pumped by a classical field with a phase (f)i, {i — 1,2,3). The 
setup could be realized in a variety of cavity-quantum-electrodynamics (cavity-QED) technologies including photonic 
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FIG. 1: (color online). Schematic representation of three interacting cavity-atom systems (^i, S2, S3) based on a possible 
implementation using photonic crystals (for illustration purposes only): the connecting wave guides carrying the driving 
classical fields with phases 01, 02, 4>3 Stxe replaced by fibers or stripline microresonators for different implementations 0, 
The three wave guides and three driving fields are labeled with the same indices to the phases (pi, 4>2, 4>3- 

crystals, circuit QED, toroidal cavities connected through fibers, Fabry-Perot cavities and coupled defect cavities 
interacting with quantum dots 0, @1- Light from the connecting waveguides/fibers can directly couple to the photonic 
modes of the atom-cavity systems through tunnehng or evanescent couphng. In each atom-cavity site we assume the 
interaction and the corresponding nonlinearity to be strong enough with at most one excited polaritonQ. 
The Hamiltonian describing the system is 

Hq =Ha,0 + Hp,o + Hjfi, (1) 

3 3 
Hafi ^"^UJcialat, Hp^ ^^^LOp^^P} P^, (2) 

1=1 i=l 



i=l i=l 

where ffa.o a-nd i/p.o are the free Hamiltonians of the wave guides and cavities, with a|, a; the field operators of 
the single- mode wave guides and tOc^i {ujp,i) the frequencies of zth waveguide mode (the polariton in ith cavity). P^ 
(Pi) the operators describing the creation (annihilation) of a mixed atom-photon excitation (polariton) at the zth 
cavity-atom system (P4 = Pi). The first summation in Hjo describes coupHngs between cavities and wave guides, 
with Ji the coupling strength between the photon mode in the zth waveguide and the adjacent two polaritons. The 
second summation in Hj o describes the classical driving of the wave guides, where ai is proportional to the amplitude 
of the ith driving field with (pi its phase and uJd the frequency of the driving fields. 

It can be seen that the Hamiltonian Hq in Eq. ([T]) is explicitly time-dependant. To remove the time dependence, 
we make the following transformation 

r)TI 

H^UlH^U^-iUl^, (4) 

where Ui = e~**"''(^i=i('^i'^»+^i^^'^. After a straightforward calculation, we obtain 

H ^ Ha + Hp + Hj, (5) 
3 3 



Ha = ^(wc,i - uJd)alai, Hp = ^(t^p,i - uJd)PiPi, (6) 

i=l i=l 
3 3 



1=1 



The density matrix p{t) of the system associated with H is related, in the following way, to the density matrix po(0 
of the system associated with Hq. 

p{t) = uipom^. (8) 



We say that the new Hamiltonian H is written in the rotating frame of the driving lasers. 
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III. THE DYNAMICS OF THE SYSTEM 

In this section, we will derive the dynamical equation for the system. 

The polaritons and waveguide modes in our system described in the last section are assumed to decay with rates 7 
and K respectively. The master equation for the whole system density operator R is: 

^ = LaR + LpR + LjR, (9) 

LaR = -i[Ha,R]+L'^R, (10) 
LpR = ~i[Hp , R] + L'pR, (11) 
LjR^-i[Hj,Rl (12) 

where Ha, Hp and Hj are given by Eqs. ([B]), and (O respectively, and 

3 2 
L'aR = I ^(2a,i?aJ - aja.i? - Raja^), L'^R = | ^(2a,i?aJ - a\c7,R - Raja,). (13) 

i=l i=l 

We use the projection operator method in Ref. Q. To this end, we define the projector PR = rss <8) tra^, 02,03-??, 
where r^g satisfying LaVss = is the equilibrium state of the three wave guides, which is close to the vacuum state 
|000)(000| when weak driving for the wave guides is assumed i.e. Ui < Ji ^ k, [i ^ 1, 2, 3). The orthogonal complement 
of P is Q = 1 — P. The operators P and Q have the properties that [13] 

PLp = LpP, (14) 
PLa = LaP = 0, (15) 
PLjP = 0. (16) 

Applying P and Q respectively to Eq. ^ and using the properties ([H]). (ITSl) and p^ . we get 

P^ = PLpPR{t) + PLjQR{t), (17) 
= QiLa + Lp + Lj)QR(t) + QLjPRit). (18) 



Formally integrate ([TSl) to get 



QRit)= / Qe^^-+^^+^'^'~'-''^LjPR{t')dt', (19) 



which is then replaced into Eq. ((T5)) . For the case Ji <C k, (/ = 1, 2,3) we only keep the second order in Ji . By 
tracing out ai, 02 and 03, we obtain 

^ = ^i[Hp,p{t)]+L'^p{t) 

POO 

+ / di'tr,,,,,,,3[Lje(^»+^'')*'Lje-^-*'(r,,®p)] (20) 







Substituting Pa, Lp and Pj with expressions p^ . pT|) and ([T^. we get 



— = - i[Pg5,p] + ^(r,_iz^_i + rj2;0Pi^,p 



Er^(^M+iP + ^+i,.P)' (21) 



j=l 



with i/gff = ^(^p,, - Wd)i^^P. + J2 ^^y^iP^P^+i + P}P^+l) + I](r»y.i^^P.+i +r,x,(/^^ + i^^i)) + h.c. , where /i.e. 



denotes the Hermitian conjugation of its previous summation. The first two summations in iJgfj cancel with each other 
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with a proper choice of Wp,, . F^^ip) = 2P,pP^. - Pjp- pP} Pj , T., = 2jf k/ {k^ +AA^), x, = -a,e''^'(2A, + iK)/(JjK), 

Hi = -2Ai/K, Ai = - i^^p,i + ujp,i+i)/2, Wp_4 = LVp^i, Zi = I + 7/(4rj), Tq = Ta and zq = z^- It can be seen 
from Eq. (|2T|) that the couphngs and detunings between the wave guide and its adjacent two polaritons induce an 
effective interaction between them given by TiUi (see i^eff)- The driving on the wave guides is equivalently transferred 
to the driving on the polaritons {TiXi in H^s), which decay with rates Ti_iZi_i + TiZi = Ti_i + + 7. Since is 
related to k, the polaritons effectively have two different channels for the decay. They can either decay directly to 
the surrounding with 7 and they can also dissipate energy via the coupling Ji_i or Ji {Jq = J3) to the adjacent two 
leaky wave guides (who also decay by k). We notice that the second channel also mixes the polaritons' operators, as 
seen in the second line of Eq. ()2ip . This mixing is actually one of the main reasons for entanglement creation among 
the polaritons. Note that the other two contributing factors are the interactions among polaritons and the driving on 
them. 



IV. COHERENT CONTROL OF THE STEADY-STATE ENTANGLEMENT 

We now derive the steady state pss by requiring that -J-11 = in Eq. (I2ip . This is done numerically due to 

. . . 

the large number of coupled equations involved. For a three-polariton density matrix, we trace out the polaritonic 
degree of freedom of cavity 1 and calculate the polaritonic entanglement of formation between cavity 2 and 3 using 
the concurrence as a measure [ll|. The concurrence C{pss) is effectively a function of the parameters Xi , yi and 
Zi. We perform a numerical optimization of C{pss) by varying these parameters and find that C{pss) is larger when 
F2 ^ Fi = F3 and 0:3 — —xi, i.e. the first and third driving fields have equal intensity but opposite phases. We also 
note here that the relation F2 <C Fi = F3 indicates that the coupling between the two cavities in question is much 
weaker than the coupling between each one of the cavities and the third cavity. Also the state of the polariton in cavity 
1 for the maximum entanglement point is found to be almost a pure state at ground energy level and therefore almost 
uncorrelated to the polaritons in cavity 2 and 3. Thus, the total density matrix p « | ground) ( ground | (8)/ci2, 3- Although 
this result initially looks counter-intuitive, it can be explained as follows: the maximum entanglement between the 
two parties, i.e. cavities 2 and 3, in a three- party system, is attained when the state of the third party, i.e. cavity 1, 
nearly factorizes in the combined three-party state. The fact that this is happening for relatively strong couplings of 
J12 = Ji and Ji3 = J3 compared to J23 = J2 is reminiscent of the behavior of a coherent process taking place. It is 
interesting to observe an analogy here with the case of coherently superposing two initially uncoupled ground states 
in a A- type quantum system through an excited state using two classical fields to mediate the interaction [l^ [l3| . 

In figure [2j we compare our setup for entanglement control of three-coupled-cavity system with the coherent 
population trapping in a three-level atom. For the latter, if the two driving fields have opposite phases and the 
atom's initial state is (|2) -I- |3))/-\/2, there will be no population in the excited state |1) and the atom will remain in a 
superposition of the states |2) and |3). Note that the states |2) and |3) are not coupled in this case. The superposition 



of them is established by a quantum interference in the state |1) 12|. It appears that the quantum correlation in 
our setup is somewhat "trapped" in the cavity 2 and cavity 3 if the driving fields 1 and 3 have opposite phases (The 
cavity 2 and cavity 3 are almost uncoupled. In this case, it is numerically verified that the driving field between them 
has almost no infiuence on their steady-state entanglement). 




|3) 



FIG. 2; (color online). The coherent trapping of a A-type three-level atom driven by two classical fields on resonance, where 
oji and ll)2 are the frequencies of the two driving fields . If the states |2) and |3) are degenerate, one could use two laser fields 
with different polarizations to distinguish the two driving paths (|2) to and |3) to |1)). 

The observation in the above paragraph is further justified by noticing that C [pss ) is varied with the phases of the 
first and third driving fields. In Fig. [3] we plot C{pss) as a function of the phases of driving fields with zi = z^ = l.OI 
and Z2 = 11. When the phase difference is (j)! ~ 4'3 = (2A; -|- l)7r {k is an integer), we get a maximum of 0.417. For 
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FIG. 3: (color online). The concurrence between the polaritons in cavity 2 and cavity 3 as a function of 4>i and 03. xi — 1.67e"^i , 
X3 = 1.67e"^^. When 01 — 03 = {2k + 1)tv {k is an integer), the concurrence reaches a maximum of 0.417. The upper left figure 
is the sectional view at 03 = 0. 




FIG. 4: (color online). Schematic diagram of the two coupled defect cavities in which there are three wave guides carrying 
the three respective classical laser fields. Note that each waveguide carrying classical fields can also be replaced by fibers or 
stripline microresonators for different implementation technologies jTi, i8j. 



general phase relations, an oscillatory behavior characteristic of the expected coherent effect takes place. There is 
a corresponding oscillatory behavior for the A-type three-level atom: the summation of the modulus square of the 
amplitudes in the states |2) and |3) is a periodic function of the phase difference between the two driving fields and 
takes a maximum when their phases are opposite • 



V. AN ALTERNATIVE SETUP: TWO COUPLED CAVITIES WITH THREE 

DRIVING FIELDS 

In Section HVl we find that when the entanglement between the two of the three cavities reaches a maximum value, 
the third cavity nearly decouples from the two cavities. It therefore seems that the third cavity plays absolutely no 
role in the establishment of the entanglement between the other two cavities. To check if this argument is correct 
and identify the role of the third cavity in the entanglement generation and control, we remove the third cavity and 
investigate the entanglement of the remaining two cavities. This new setup is shown in Fig. |3J where there are 
three wave guides coupled to two cavity-atom systems and these three wave guides are driven by three classical fields 
respectively. We analyze the polaritonic entanglement between cavity 2 and 3 (relabeled as Si and 5*2 in Fig. S]). 

The Hamiltonian and the derivation of the effective master equation are similar to those for the three-cavity setup 
in Section |TT] and IIIII We therefore omit the detailed derivation steps and provide only the final effective master 



FIG. 5: (color online). The concurrence between two cavities -FigUJ- as a function of (^i and (^3. X2 = 1/2 = 0, xi — 5e'^''- ,xs = 
5e"^^, Fi = Fa = 1.316 x 10* and F2 = 10^". When 4>i ~ 4>'i ~ (2fc + l)7r {k is an integer), the concurrence reaches a maximum 
of 0.470. The upper left figure is the sectional view at (^3 = 0. 



-Qie*'^'(2Ai+zK)/(JjK), y2 = -2A2/K, Ai = Wc,i-Wp,i, A2 = a;c,i-(wp,i+a;p,3)/2, A3 = a;c,i-Wp,3, Zi = l+7/(4rj). 

The optimization of this entanglement gives similar values of the parameters like the ones used above except that 
the values for are reversed, i.e. r2 S> Fi = F3; however, the concurrence reaches a maximum of 0.47. Again the 
dependence 4>i — (f)^ = {2k + l)7r (fc is an integer) is apparent (see Fig. [5]). However, if we compare the insets in 
Fig. [3] and Fig. [5] for the cross-sectional plots of the concurrence for (f)j, — 0, we see that the plot in Fig. [3] has a 
narrower peak whereas the plot in Fig. [S]is broader. This implies that the maximum concurrence for configuration in 
Fig. [4] is substantially more stable against variation in the phases and (p^ than that in Fig. [1] However, when the 
dissipation (parametrized by 7 in Zi)) increases, the entanglement in the latter configuration decreases more slowly 
than the former one. This can be numerically verified. Thus we conclude that cavity 1 in Fig. [l]not only mediates 
coherently between cavities 2 and 3, but it also stabilizes the amount of entanglement between the two cavities. 

There are many other configurations for the coupled-cavity setup. For instance, one could consider an extension of 
the setup in Ref. [H to three defect cavities, as shown in Fig. |6l However, numerical optimization for this extension 
and many others does not seem to increase the polaritonic entanglement between any two cavities. Therefore, the 
setups in Fig. [T] and S] appear to be optimal ones for two-polariton entanglement. 



equation. 



dp 



- ^I^cff'/'] 

+ {T2Z2 + ri)FfiP + (F2Z2 -t- F3)F2^2P 



dt 



(22) 



with H'^Q = (F2i/2-Pi-P2 + ri+i2^j+i)-P/) + h.c. , whcrc h.c. denotes the Hermitian conjugation of its 

previous summation. F^j{p) is defined in Section [Hi] as 2PipP^ — P-Pjp — pPjPj , Tj = 2jfK/(K^ + 4A^), Xi — 




FIG. 6: (color online). Three defect cavities coupled to one wave guide. 
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VI. THERMALIZATION OF THE COUPLED-CAVITY SYSTEM 



In this section, we consider the thermalization of the lossy driven atom-cavity system. For simplicity, we consider 
a simpler system which involves two defect cavities coupled to a driven wave guide, as shown in Fig. [T] 




FIG. 7: (color online). Two defect cavities coupled to one wave guide. 

This system was studied in Ref. [sj, where the reservoir temperature is set to be zero and an analytical solution 
was obtained (see Eq. (23)-(29) therein). For finite temperature, the master equation needs to be modified, i.e. Eq. 
(12) and (13) of Ref. Q are replaced by 



L'^R =K{nc + l)(2ai?at - a^aR - Ra^a) + Knc{2a) Ra - aa) R - Raa)), 

2 

L'pR^^-f{np + l){2a^Ral - aja^R - paja^) + jnp{2ajRa - aajR - Ra^aj), 



(23) 
(24) 



i=l 



is the mean photon number at the reservoir temperature Tr and the cavity frequency 
is the mean photon number at the reservoir temperature Tr and the polaritonic 



where Uc — 
LOcav Similarly Up 
frequency ujpoi- 

The effective master equation for the two polaritons can be obtained using the same method in Ref. [3|. For 
Tfl ^ hupoi/ks, the temperature terms in Eq. (j24p are preserved in the final effective master equation i.e. Eq. (20) 
of Ref. Q. The steady state p"^ is obtained by requiring = 0. To characterize the degree of thermalization of 
the steady state, we calculate the distance between the steady state and a thermal state, using the following distance 
measure 



lysta 
0: 



d{p' 



„th\ 



1 



th\ 



(25) 



The trace distance d{p^'^ , p ) provides a useful measure to distinguish the steady state p"" from the thermal state 
p*'' through quantum measurements ^5j. Therefore, if d{p^'' , p*'*) increases with system parameters we say that the 
system is farther away from thermalization. Also, the thermal state p*'' is chosen to be p*'* = eyi\)[—fiujpoi{(T\(Ji + 
^2'^2) /kBTn] up to a normalization factor tr(p*''). 

as a function of x 



Fig. |8(a) shows the distance d{p'^'' , p*'*) 



and Tr, where x is a parameter defined in Ref. Q 



(below Eq. (22)) and it is proportional to the strength of the driving field 
(see Ref. j3|) 



The relevant parameters y — 15, z = 1.01 



The unit of is hupoi/kB- It is seen in Fig. 8(a) that the steady state is close to the thermal state 
if there is no driving field, and for stronger driving field the steady state is farther away from thermalization. This 
is reasonable from a physical perspective as the driving field generally induces coherence (i.e. non-zero off-diagonal 
elements in the polaritonic density matrix) for the polaritons while the thermal state is diagonal. In addition, it seems 
that d{p'"', p*^) does not depend on the reservoir temperature. This may be because Tu <C hujpoi/kB so that the 
effect of the thermal agitation is rather small. The effect should certainly manifest itself for larger T/j. However this 
regime is beyond the approximation for the derivation of the effective master equation (Tr ^ hujpoi/kB) and it is in 
general not easily solvable even with numerical calculations. 

Comparing Fig. |8(a) for a fixed Tr with the first plot of Fig. 2 {y = 15) in Ref. i3:|, one finds that they are not 
consistent, especially for large x, for which d{p 
However, if one takes the derivative of d{p''^ , p 



p"'") is very large while the polaritonic entanglement is negligible. 
'*'') with respect to x, then a relationship appears. Fig. |8(b)| shows 
\dd{p''^ , p*'^)/dx\ as a function of x and Tr. It can be seen that there are two peaks for a fixed temperature. 
This is similar to the first plot of Fig. 2 in Ref. [1]. Also the two plots are consistent for large x. Therefore, it 
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FIG. 8: (color online). The distance d(p°'' , p***) for (a) and the derivative |9d(p''^, p*'')/92;j for (b) as functions of a; (proportional 
to the driving strength) and Tr. 



may be concluded that the change rate of the thermalization with respect to the driving strength (rather than the 
thermalization itself) is related to the polaritonic entanglement. Physically, for a increase/decrease of the driving 
strength i.e. more/less coherent energy is injected into the system, a more rapid change of the thermal property (or 
the degree of thermalization) of the system indicates that a stronger correlation (entanglement) is established. The 
coherent energy refers to fact that the driving field induces off-diagonal elements in the polaritonic density matrix as 
mentioned previously. One could conjecture that a more rapid change of the degree of thermalization of the system 
may indicate that the interaction between the two polaritons are stronger which leads to a stronger entanglement 
between them. 



VII. CONCLUSION 

In this paper, we show that long-distance steady state entanglement in a lossy network of driven light-matter 
systems can be coherently controlled through the tuning of the phase difference between the driving fields. The 
role of driving phase field in engineering interaction and entanglement in coupled atom-cavities was also discussed 
in Ref. [l6|. Here, it is found that in a closed network of three-cavity-atom systems the maximum of entanglement 
for any pair is achieved even when their corresponding direct coupling is much smaller than their couplings to the 
third party. This effect is reminiscent of coherent effects found in quantum optics that coherent population transfers 
between otherwise uncoupled levels through a third level using two classical coherent fields. An alternative geometry: 
two-coupled cavities with three driving fields is discussed. For finite temperature, we analyze the thermalization of 
the two defect cavities coupled to one driven wave guide. It is found that the change rate of the thermalization of 
the system with respect to the driving strength (rather than the thermalization itself) can indicate the degree of the 
polaritonic correlation (entanglement). 
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